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WE CLAIM: 

1. A method of operating a plurality N of seismic vibrators 
simultaneously witli continuous sweeps, and separating the seismic response 
5 for each vibrator, said method comprising the steps of. 

(a) loading each vibrator vi^h a unique continuous sweep 
signal consisting of M ^ N segments, the /* segment being of the same 
duration for each vibrator, / « 1, 2, . . . , M; 

(b) activating all vibrators and using at least one detector to 
10 detect and record the combined seismic response signals from all vibrators; 

(c) selecting and recording a signature for each vibrator 
indicative of the motion of that vibrator; 

(d) parsing the vibrator motion record for each vibrator into M 
shorter records, each shorter record coinciding in time with a sweep segment, 

15 and then padding the end of each shorter record suffidently to extend its 
duration by substantially one listening time; 

(e) forming an M x N matrix s whose element Sff(t) is the 
padded shorter vibrator motion record as a function of time / for the vibrator 
and 7^ sweep segment; 

20 (f) parsing the seismic data record from step (b) into M 

shorter records, each shorter record coinciding In time with a padded shorter 
record of vibrator motion from step (d); 

(g) fomning a vector d of length M whose element <// is the 
shorter data record from the preceding step; 
25 (h) solving for Ej(/) the following system of M linear equations 

in iV unknowns 

SE=D 

where S(f(f) Is the Fourier transfonm to the frequency (f) domain of S(/(t) and 
2>i(/) is the Fourier transfonn of rf^O. where / = 1, 2, . . . M and> = 1, 2, . . . 
30 N; and 

(1) Inverse Fourier transforming the E/f) to yield e/O- 
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2. The, method of claim 1, wherein each sweep segment is 
seiected from one of the following sweep-design categories: (a) linear, (b) 
nonlinear, and (c) pseudo-random. 

3. The method of claim 1, wherein all of the N unique continuous 
5 sweeps are identical except for the phase of their segments. 

4. The method of claim 3, wherein all N segments are identical 
except for phase, and the phase differences for the iV^ sweeps are determined 
by the following steps: (a) constructing a reference sweep by starting with a 
preselected reference segment, then advancing the segment 360/Af degrees 

10 in phase to make the second segment, then advancing the phase 360/M more 
degrees to malce the third segment, and so on to generate a sweep of M 
segments: (b) constructing a first sweep by advancing the phase of the first 
segment of the reference sweep by 90 degrees; (c) constructing a second 
sweep by advancing the phase of the second segment of the reference sweep 

IS by 90 degrees; (d) and so on until JST sweeps are constructed. 

5. The method of claim 1, vdiereln each unique continuous sweep 
has a duration in time sufficiently long to collect all seismic data desired 
before relocating the vibrators. 

6. The method of claim 1, wherein the vibrator signature record for 
20 each vibrator is a weighted sum or ground force record of the motion of that 

vibrator. 

7. The method of claim 1 , wherein M=N and the system of linear 
equations SE = D is solved by matrix methods comprising the steps of 
deriving a separation and inversion filter (SI}''* by inverting the matilx S. then 

25 perfomiing the matrix multiplication (S )~^D . 

8. The method of claim 1 , wherein the system of linear equations 
SE-D\s solved by matrix methods and the method of least squares 
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comprising the steps of deriving a separation and inversion filter of the form 
F = (s'sy^S* , then perfonning the matrix muitiplication FD . 

9. The method of claim 1, wherein each segment has a duration 
that is at least as long as the seismic wave travel time down to and back up 
S from the deepest reflector of interest. 



